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We investigate in a general form entanglement of biphoton qutrits and ququarts, i.e. states formed 
in the processes of collinear and, correspondingly, degenerate and non-degenerate Spontaneous 
Parametric Down-Conversion. Indistinguishability of photons and, for ququarts, joint presence of 
the frequency and and polarization entanglement are fully taken into account. In the case of qutrits 
the most general 3-parametric families of maximally entangled and non-entangled states are found, 
and anti-correlation of the degree of entanglement and polarization is shown to occur and to be 
characterized by a rather simple formula. Biphoton ququarts are shown to be two-qudits with the 
single-photon Hilbert space dimensionality d = 4, which differs them significantly from the often 
used two-qubit model (d — 2). New expressions for entanglement quantifiers of biphoton ququarts 
are derived and discussed. Rather simple procedures for a direct measurement of the degree of 
entanglement are described for both qutrits and ququarts. 



INTRODUCTION 



Optical qutrits and ququarts are promising objects 
of the modern quantum information and quantum 
cryptography 0-0]. Formally, qutrits and ququarts are 
defined as superpositions of, correspondingly, three and 
four basis states. In practice, most often, the basis 
states for qutrits and ququarts are formed by biphoton 
states arising in the processes of Spontaneous Parametric 
Down-Conversion (SPDC). For qutrits it is sufficient 
to use the collinear degenerate SPDC processes, i.e., 
such processes in which wave vectors of two photons 
in a SPDC pair are strictly parallel to each other and 
frequencies are also given and equal to each other. For 
constructing ququarts, one has to use either the non- 
collinear frequency-degenerate or collinear but frequency- 
non-degenerate SPDC processes, i.e., processes in which 
either directions of wave vectors or frequencies of two 
photons in SPDC pairs differ from each other. Below we 
will refer both of these possibilities as "non-degenerate" . 
In theory, from the very beginning Q, biphoton qutrits 
were considered as arbitrary superpositions of three Fock 
states, corresponding to three possibilities of distributing 
two indistinguishable photons of an SPDC pair in two 
polarization modes, horizontal (H) and vertical (V) ones. 
Each of these basis states is a direct product of two 
one-qubit single-photon states, and biphoton qutrits are 
two-qubit states. There were many works devoted to 
biphoton qutrits (fj|-tl5|. Most of them study polarization 
properties of qutrits and much less entanglement. Giving 
a comprehensive picture of entanglement of qutrits 
is one of the goals of this paper. In particular, we 
will find and describe general families of maximally 
entangled and non-entangled qutrits, their entanglement 
quantifiers such as the Schmidt entanglement parameter, 



concurrence, and the subsystem entropy, relations 
between entanglement and polarization of qutrits, etc. 

An important point to be discussed is the role of 
the symmetry of biphoton wave functions with respect 
to permutation of particles' variables. Existence of 
entanglement related to symmetry was realized by 
many authors rather long ago, [lj|-[2(j] and so on. 
But a general attitude to the entanglement related 
to symmetry is rather sceptical, up to the opinion 
that this type of entanglement is unimportant and 
can be forgotten. In such approach all basis states 
of qutrits and ququarts would be non-entangled and 
the only reason for entanglement would be related to 
various choice of coefficients in superpositions of basis 
states. This is a configurational entanglement. But in 
reality, owing to symmetry, basis states of qutrits and 
ququarts are entangled (at least one basis states in 
the case of qutrits and all basis states of ququarts). 
This is a fundamental, unchangeable entanglement. 
In superpositions of basis states the symmetry-related 
and configuration entanglement exist together and each 
of them is so strongly built into a general picture 
of entanglement of qutrits and ququarts, and into 
their entanglement quantifiers, that it's impossible to 
separate these two types of entanglement and to take off 
the symmetry-related one without hurting significantly 
the general picture and results. Thus, all kinds of 
entanglement have to be taken into account together and 
none of them can be "forgotten" . 

In the case of ququarts, the symmetry-related 
entanglement arises from both polarization and 
frequency (or angular) degrees of freedom. This makes 
the traditional two-qubit model of biphoton ququarts 
invalid. As we show, biphoton ququarts are two-qudit 
states with the dimensionality of the one-photon Hilbert 
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space d = 4 and dimensionality of the two-photon 
Hilbert space D = d 2 = 16. This makes ququarts 
significantly different from qutrits (where d = 2 and 
D = d 2 = 4) , and a series of new results on entanglement 
of ququarts is derived and discussed in Section 7. This 
new understanding of the physics of ququarts raises 
a question about changes in applications of ququarts 
analyzed earlier in the frame of the two-qubit model 
when the latter is substituted by the qudit picture. We 
hope to return to such analysis elsewhere. 

In addition to the above-mentioned entanglement 
quantifiers, we use widely the Schmidt decomposition 
of biphoton wave functions [U |22] . As known (23[, 
for pure biphoton states with the dimensionality of the 
one-photon Hilbert space d the maximal amount of 
terms in the Schmidt decomposition equals d, and such 
states are unseparable. Only if all but one coefficients 
in the Schmidt decomposition are equal zero, and 
the remaining exceptional coefficient equals unit, the 
Schmidt decomposition is reduced to a single product of 
Schmidt modes, and such state is separable. We consider 
this criterion as the ultimate indication of whether states 
are separable or not. 



2. STATE VECTORS AND WAVE FUNCTIONS 
OF BIPHOTON QUTRITS 

In the form of state vectors, purely polarization 
biphoton states (qutrits), are given by a superposition 



bosons the corresponding formula has the form [24J (in 
slightly modified notations) 

*(aci, x 2 ,...,x n ) = (xi, X2,...,x n \m, n 2 ,..., n k) 

= , , ] f=? y^-P(fli(si)gi(s2)--gi(s ni ) 

xg 2 (x ni+1 )g2(x ni+ 2)--g2 
xg k (x ni+n2+ ... +nk _ 1+1 )...g k (x„)j , (2.4) 

where x\, x 2 ,---, x n are dynamical variables of identical 
boson particles, gj(xi) are single-particle wave functions 
(of j-th modes and i-th variables), P indicates all possible 
transpositions of variables Xi in wave functions gj(xi), 
m, 112, rik are numbers of particles in modes, n is the 
total amount of particles in all modes, k is the total 
amount of modes; for empty modes the corresponding 
single-particle wave functions have to be dropped. 

In the case of qutrits we have two modes (j = H or V) 
and two particles, n = k = 2. The polarization variables 
of two photons can be denoted as o\ and <J2 ■ In terms of 
wave functions, the single-photon wave functions gj(xi) 
are given by the Kronecker symbols. Thus, the qutrit 
basis wave functions corresponding to the basis state 
vectors in Eq. (|2.ip can be written as 



*ffff(CTl, CT 2 ) = (ai,a 2 \2 H ) = S^HSa^H, 

^hv (01,02) = (ai,a 2 \lH, lv) 



I*) = Ci\2h) + C 2 \1h, lv) + C 3 |2 T 



(2.1) 
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(2.6) 
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where the basis state vectors are given by 
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|0>, 



V2 

(2.2) 

|0) is the vacuum state vector, a^ H and Oy are the creation 
operators of photons in the modes with horizontal 
and vertical polarizations (with given equal frequencies 
and given identical propagation directions). are 
arbitrary complex constants Ci.2,3 = |Ci j 2,3|e l¥ ' 1 ' 2 ' 3 , 
obeying the normalization condition 



\Ci\ 



\C 2 \ 2 



\c 3 \ 2 



(2.3) 



Actually, as the total phase of the state vector (|2.ip or 
wave function (see below) does not affect any measurable 
characteristics of qutrits, one of the phases f 103 , or a 
linear combination of phases, can be taken equal zero, 
and, hence, the general form of the qutrit state vector 
(|2.ip is characterized by four independent constants (e.g., 
|Ci|, \C 3 \, (p\, and ip 3 with tp 2 = 0). 

As qutrit (|2 . 1 [) is a two-photon state, its polarization 
wave function depends on two variables. A general rule of 
obtaining multipartite wave functions from state vectors 
is known pretty well in the quantum-field theory, and for 



The same basis wave functions can be written 
equivalently in the form of 2-row columns, which is more 
convenient for calculation of matrices 
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*w = — 



■ 



(2.8) 
(2.9) 




(2.10) 



where the upper and lower rows in two-row columns 
correspond to the horizontal and vertical polarizations 
and the indices 1 and 2 numerate indistinguishable 
photons. 

In a general form the qutrit wave function 
corresponding to the state vector (|2.1[) is given by 



vv, 



(2.11) 
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where ^hh, ^hv, an d ^vv can be taken either in the 
form (pj]) - (|27Hl or CT- d^TUl) . 

Alternatively, the same general qutrit wave function 
p. lip can be presented in the form of an expansion in a 
series of Bell states 



where 



C!±C 3 

V2 



(2.12) 



(2.13) 



= * 2 (j2~TU|) and 




(2.14) 



^ and \I> + are the wave functions describing three Bell 
states. The fourth Bell state, 



1 

V2 



, (2.15) 



does not and cannot arise in the expansion (|2.12p because 
(|2.15p is antisymmetric with respect to the variable 
transposition 1^2, whereas all biphoton wave functions 
have to be symmetric. Nevertheless, in principle, the 
antisymmetric Bell state (|2.15p can be included into the 
expansion (|2.12p but obligatory with the zero coefficient: 



(2.16) 



This obligatory zero coefficient in front of \I r_ or missing 
forth antisymmetric Bell state in the expansion (|2.12|) 
is related to restrictions imposed by the symmetry 
requirements for two-bozon states. Thus, even existence 
of biphoton qutrits as superpositions of only three basis 
wave functions occurs exclusively owing to the symmetry 
restrictions eliminating the forth (antisymmetric) basis 
Bell state. 



3. DENSITY MATRICES 

The first step for finding the degree of entanglement 
is related to a transition from the wave function ^ to 
the density matrix of the same pure biphoton state p = 
"J"!^. The full density matrix of the qutrit (|2.1ip can be 
presented in following two forms 
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(3.1) 



and 



|Ci| 2 
75 Cl * C2 

1 /nf* /-» 



72 ClC2 73 ClC i 



±|C 2 | 2 i|C 2 | 



ClC; 



C3C2 



C*3*C 2 



(3.2) 



|C 3 | 2 / 



The next step is the reduction of the density matrix with 
respect to one of the photon variables, e.g., of the photon 
2. Mathematically this means taking traces of all matrices 
with the subscript 2 in Eq. (|3.ip . which gives 



Tr 2 p 



c : 1 

C\C2 + C2C3 



|C 2 | 2 CiC 2 * 



V2 



1 a 



V2 

12 , l 6 2 



C2C|\ 



O' 2 



(3.3) 



It may be interesting to analyze a relation between the 
4x4 density matrix (|3.2p and the 3x3 coherence matrix 
introduced by Klyshko in 1997 (25|. The density matrix 
p (|3.2p is written in a natural two-photon basis 



(3.4) 



The question is how it can be transformed to the basis of 
states ^hh, ^>vVi ^ hv plus the empty antisymmetric 
states <f>~ (|2.15p ? Evidently, the transformation 





(3.5) 
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is provided by the matrix 



U = 



(l 









1 

v/2 
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y/2 
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Now. transformed to the basis {^hh , ^vv-, 
the density matrix p (|3.2p takes the form 



HV 



Ptransf = U pU 





Vcrc 3 



IC 2 I ; 



c 3 c. 



dC* z \ 




o |c 3 |V 



(3.6) 



(3.7) 



A part of this matrix with nonzero rows and columns 
coincides with the 3x3 coherence matrix [H, [25[ 



Pcoh — 



(\Ci 

C* C3 I C 2 

\ctc 3 C Z G* 2 



2 c 3 *c 2 



(3.8) 



|C a 



Though the coherence matrix p co ft (|3.8p is widely used 
and analyzed in literature, both p co h and ptransf are 
hardly appropriate for reduction over one of the photon 
variables (e.g., 2) and for finding correctly the reduced 
density matrix p r (|3.3p because the variables 1 and 2 are 
mixed not only in the matrix ptransf itself but also in the 
transformed basis of Eq. (|3.5[) . 



4. DEGREE OF ENTANGLEMENT 

As known (2~H [22j], the trace of the squared 
reduced density matrix p 2 , determines purity of the 
reduced state coinciding with the inverse value of the 
Schmidt entanglement parameter it -1 . The result of its 
calculation for the reduced density matrix of Eq. (|3.3[) is 
given by 



K~ 



Tr{pl) 



dl 



ic 2 i 2 



c. 



\Co 



+\ctc 2 + c* 2 c s \ 



(4.1) 



With the normalization condition (|2.3[) taken into 
account, Eq. (|4.1|1 can be reduced to a much simpler form 



K 



|2CiC 3 



(4.2) 



It's known also |23|, that in the case of bipartite states 
with the dimensionality of the one-particle Hilbert space 
d = 2, there is a simple algebraic relation between 
the Wootters' concurrence C (2(| and the Schmidt 
entanglement parameter K, owing to which 



C 




\2CiC z - C\ 



(4.3) 



At last, in terms of the constants C± l[2.13|) . Eqs. (H^J) 
and f|4. 3[) take the form 



K = 



2 - \Cl - Ci 



-gtp, C=\C\-C^-Cl\. (4.4) 



Note that the expressions for the concurrence C [Eq. 
(|4.3|) and the last formula of Eq. (|4.4[) ] can be derived 
also directly from the original Wootters' definition (2(|. 
Indeed, for a pure bipartite state with the dimensionality 
of the one-particle Hilbert space d = 2, the concurrence 
is defined in Ref. [2g| as 



c = k*!** 



(4.5) 



where \& is the function or state vector arising from W 
after the "spin-flip" operation 



* = K)iK) 2 *. 



(4.6) 



and a y is the Pauli matrix, a y 



J). For qutrits, ^ 
is given by Eqs. (|2.1ip or (|2.12p . The rules of the "spin- 
flip" transformation for one-photon wave functions are 
( ) — » — t (j) and (i)^*(o) - From here we easily find 
the spin-flip transforms of the qutrit basis wave functions 

$+ = = (4.7) 

and of the general-form qutrit wave function (|2.11l) 



= -(7+*+ + C_*_ + C 2 ^hv- 



(4.8) 



Substitution of these expressions into Eq. (|4.5|) gives 



C 2 = 



-2C\Cl+C* 2 



c 2 * 



(4.9) 



in a complete agreement with Eqs. (|4.3[) and (|T.T|) . 

In a special case of real constants Ci,2,3 an d 
C±, owing to normalization (|2.3p . the Schmidt 
entanglement parameter and concurrence (|4.4p appear to 
be determined by the only real parameter C+ 



1 + 4C| - ac: 



\2Ci 



1 



(4.10) 



The functions K rea i(C+) and C rea ;(C+) are shown in 
Fig. [T] together with the subsystem entropy found in the 
following section. This picture demonstrates clearly that 
qutrits are non-entangled (K rea i = 1 and C rea i = 0) if 
C + = ±l/\/2 and, hence, Ct + Cjj V = \. Consequently, 
the family of non-entangled qutrit wave functions with 
real coefficients is given by 



^NEreal = ( X 



COS 



HH 



V2 



sm 



*w (4.11) 
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-1 -1/V2 -0.5 



C + 0.5 1/V2 1 



Phc. 1: The Schmidt entanglement parameter K rea i 14,20 . 
concurrence C rea ; l |4,3|l and the von Neumann subsystem entropy 
S rrea i 15,71 of the qutrit 12. lit , 12.121 1 with real constants Ci,2.3, 
C± vs. C+ 



with arbitrary </>. If the constants 61,2,3, C± are complex, 
a general condition of no-entanglement is the same: 
C = 0. As seen from the general expression for C 
(|4.3j) . the concurrence depends on phases (^1,2,3 of the 
constants £1,2,3 only via the combination ipi + 993 — 2(^2 ■ 
Hence, if 932 = 2 "( < ^ 1 m the case of complex 



constants 6*1,2,3, Eq. (|4.3j) is reduced to C 



|C 2 | 2 



2|Ci||C 3 | 



i.e. this case appears to be equivalent to the 

case of real constants Ci,2,3- From here we find that the 
general three-parametric family of wave functions of non- 
entangled qutrits is given by 



*NE g eneral(<t>,Vl,<P3) = ^ e^+^tf 



HV 



cos 2 ( I) e^HH +sin 2 I | ) e^ vv (4.12) 



with arbitrary (f>, tpi, and (^3. With these wave functions 
found explicitly we can reconstruct the corresponding 
family of state vectors of non-entangled qutrits 



V&jNE general = 



£(¥>i+¥>3)„t „t 
H a V 



(|) e^og +sin 2 (0 e^} |0). (4.13) 



Qutrits are maximally entangled when C = 1 and 
A' = 2, and for the wave functions with real constants 
61.2.3, C± this occurs in two cases: C+ = ±1 and 
C+ = 0. In the first of these cases ^ = ±\&+. In the 
second case the maximally entangled wave function has 
a form of an arbitrary superposition of ^ hv and <I>_, 
^ = smcfi^HV + cos</>\I/_. As previously, this result 
can be generalized for the case of wave functions with 
complex coefficients such that ip2 = j(v>i + V3)- As a 
result we get the following three-parametric family of 



wave functions of maximally entangled qutrits 
*max(0, <pi , <^ 3 ) = sin e^+^HV 

v2 

and the corresponding family of state vectors 

-{' 



(4.14) 



COS( 



^ 2 )}|o>. 



(4.15) 



Eqs. (|4.14p and (|4. 15)) are more general than the 
expression for the maximally entangled state of Ref. 
[27j . which follows from (|4.15| at <f> = 0, |*max)U=o = 
^(e^|2ff)-e 



' 3| 2v»- 

= tt/2 Eqs. (gmj) and (j^7T3|) show 



In the case </> 
that the state ^hv (\^h,1v)) belongs to the family 
of maximally entangled states too (contrary to a met 
opinion that the state \1h> lv) is factorable). A way of 
seeing explicitly whether qutrits are factorable or not 
consists in finding their Schmidt decompositions, which 
can contain either two products of Schmidt modes or only 
one product. This analysis is carried out in the following 
section. 

But before switching to the Schmidt-mode analysis let 
us discuss briefly the problem of qutrit polarization. If we 
define the biphoton polarization vector, as suggested by 
Wang [2§ |. £ = Tr (p r (?), where a is the vector of Pauli 
matrices, from Eq. (|3 . 3[) we easily find 

f = {V2Re(CiC 2 * + C 2 C 3 *), -V^lmidCI + C 2 C 3 *), 
|C X | 2 - |C 3 | 2 }- (4.16) 

A direct comparison with the results obtained in 1999 by 
Burlakov and Klyshko Q for polarization characteristics 
of qutrits shows that the polarization vector £ (|4. 16|) 
coincides exactly with one half of the vector of Stokes 
parameters S = {Si, S2, S3} of Ref. ||, and the absolute 
value of £ coincides with the degree of polarization P 
(introduced also in Q) 



1 



S. 



= P = 



Sj 



si 



si 



(4.17) 



Note that the Stokes parameter and the degree of 
polarization were found in Ref. Q in a way, absolutely 
different from that used above for derivation of the 
polarization vector £. For this reason the coincidences 
(|4.17p are rather non-trivial. 

The derived expression for the polarization vector £ 
(|4.16j) can be compared with the general expression for 

the qutrit concurrence C (|4.3[) to show that they obey 

-.2 „ 

the relation £ = 1 — obtained by Wang in 2000 
[28[. Note, however, that for qutrits we have used the 
concurrence of Eq. (|4.3[) . C = \2C\Cz — C||, rather than 
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2\ad — bc\ obtained in Refs. |28| for the state |^) 6cd = 
a|00) + fo|01> + c|10) + d\ll) to be commented below in 
section 7. In terms of the degree of polarization P the 

relation £ = 1 — C can be rewritten as 



C 2 + P 2 



1. 



(4.18) 



Thus, the degrees of polarization and entanglement 
anti-correlate with each other: the maximally entangled 
qutrits are non-polarized, and maximally polarized states 
are non-entangled. Though, maybe, intuitively more or 
less expected, as far as we know, anticorrelation of 
polarization and entanglement has never been presented 
in a rigorous mathematical form of Eq. (|4.18|) . 



5. SCHMIDT MODES OF QUTRITS AND THE 
SUBSYSTEM ENTROPY 

Entanglement means that the biphoton wave function 
cannot be factorized whereas no-entanglement means 
that it is factorable. A transition from non-factorable to 
factorable wave functions can be reasonably explained 
in terms of Schmidt modes. Schmidt modes are 
eigenfunctions of the reduced density matrix, i.e., 
solutions of the equation p r ip = XtJj. As in the case of 
qutrits p r (|3.3p is the 2x2 matrix, it has two eigenvalues 

A±, its eigenfunctions are 2-row columns ip± = 

and the eigenvalue- eigenfunction equation has the form 



«± 
b± 



(5.1) 



The modes can be normalized |a±| 2 + \ b±\ 2 — 1 and they 
are orthogonal to each other a+a*_+b + b*_ = 0. In terms of 
A± the Schmidt entanglement parameter equals to K — 

(x 2 + + x 2 _y\ 

In accordance with the Schmidt theorem, the biphoton 
wave function can be presented as a sum of products of 
Schmidt modes (Schmidt decomposition) . In the case d = 
2 (two-qubit states or qutrits) the Schmidt decomposition 
contains only two terms 



£V*Ilk=(lW±(2), 



(5.2) 



where arguments of the Schmidt modes indicate variables 
of photons "1" and "2" . This decomposition shows that 
in a general case the wave function is nonseparable. 
Exceptions occur when one of the eigenvalues of the 
reduced density matrix, A+ or A_, becomes equal zero. 

Eigenvalues of the matrix (|3.3[) can be found rather 
easily and they can be reduced to a very simple form 
being expressed via the concurrence C 



Ah 



(5.3) 



In the case C = 1, when entanglement is maximal, 
A + = A_ = i, i.e., two products in the Schmidt 



decomposition (|5.2[) are presented with equal weights 
and, clearly, the wave function is nonseparable. 

In the case C = (no entanglement) Eq. (|5.3p gives 
A+ = 1 and A_ = 0. In this case one of the products of 
Schmidt modes (?/>_?/;-) in the Schmidt decomposition 
disappears because the coefficient \/X~Z in front of 
it vanishes. This is the reason of separability of the 
wave function in the case of no-entanglement. As for 
the remaining product in the Schmidt decomposition, 
the corresponding eigenfunction of the reduced density 
matrix ip + at C = can be found easily and has a 
reasonably simple form (|4.12|) 



'cos(0/2)e* Vl 
sin(0/2)ei^ 3 



(5.4) 



With this expression we find immediately that in the 
case C = the Schmidt decomposition (|5.2[) gives the 
following representation for the factored biphoton wave 
function 

= / cos(0/2)e^ \ / cos(^/2)e^ \ 
NE ^sin(0/2)e5^ J ^sin(0/2)e^3 J ' [ ' J 



Identity of Eqs. (|4.12p and (|5.5p can be checked directly, 
be means of substitution in Eq. (|4.12p instead of ^hh, 
^hv, ^vv their column representations (|2.8|) - (|2.10l) . 
calculation of direct products of columns in all terms of 
Eq. (|4.12jl and in Eq. ()5.5j) , and presenting ^ ne in both 
cases in the form of the following 4-row column 



^NE general (0, <Pl > ^3) 



/ cos 2 ((f>/2)e^ 1 \ 



(5.6) 



\ sin 2 (</>/2)e iV3 J 
Eq. (|5.3p can be used to find the subsystem entropy 



123 defined as 



S r = -Tr (p r log 2 pr) = - ^ A± log 2 



(5.7) 



In the case of qutrits with real coefficients the 
concurrence C rea i itself is a function of C+ (|4.10j) . With 
C rea i(C + ) substituted instead of C into Eq. (|5.3[) and 
then A±[C(C+)] substituted into Eq. (|5.7[) . we get a 
function S rrea i(C+), which is plotted in Fig. [TJ Though 
different from both C(C + ) and K(C + ), the subsystem 
entropy S rr eai(C+) has the same main features as two 
other entanglement quantifiers: S rrea i is minimal and 
equals zero at C+ = ±l/\/2 and S r rea i is maximal and 
equals 1 at C+ = and C+ = ±1. Therefore, though 
the Schmidt entanglement parameter K, concurrence C 
and the subsystem entropy S r characterize the degree of 
entanglement in different metrics, their behavior is very 
similar, which confirms all conclusions made above about 
conditions of separability and nonseparability of qutrits 
and their wave functions. 
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As mentioned above, one of the basis wave functions of 
qutrits, ^ hv (|^-6|) . is maximally entangled (K = 2, C = 
1) and, hence, unseparable. Eigenvalues of the reduced 
density matrix of this state are degenerate, A+ = A_ = i , 
and the Schmidt decomposition has the form 



^ + (l)V+(2)+^_(l)V-(2) 



with the Schmidt modes given by 



Mi 



= ( ; ) , ^ = 



-f 1 
V2 v-i 



(5.9) 



Thus, entanglement of the state ^ hv Q2.6P is proved 
here in several ways, by direct calculations of the 
entanglement quantifiers K 1 C and S ri and by showing 
that its Schmidt decomposition contains two products of 
the Schmidt modes (|5.8p . In literature the opinion about 
entanglement of the state \& hv is shared by some authors 
E3, HOI ; thou gh s ometimes the state ^hv is claimed to 
be separable [301 ■ Actually, this difference of opinions 
reflects a rather popular point of view that all basis states 
of qutrits and ququarts (see below) are separable and 
non-entangled. In a more general form, this statements 
can be reformulated as saying that all states of the type 
a ! a ]|0) generate separable wave functions in both cases 
of coinciding i = j and non-coinciding i ^ j mode 
indices. Our analysis shows that for biphoton states the 
latter is not true. In reality, for arbitrary photon variables 
(polarization, angular, frequency, or their combinations), 
if ipi and tjjj are the wave functions of the i — th and j — th 
modes and i ^ j, the state vector oja||0) generates the 
wave function of the type ^ hv lj2.6[) . 



^(1)^(2) + ^(2)^(1) 
V2 



(5.10) 



which is entangled with the degree of entanglement 
characterized by K = 2, C = S r = 1, and the Schmidt 
decomposition of the type (|5.8|) , ([52 



6. FINDING THE DEGREE OF 
ENTANGLEMENT OF QUTRITS FROM DIRECT 
POLARIZATION MEASUREMENTS 

In this section we will show that there is a rather 
simple method of measuring experimentally the degree 
of entanglement of qutrits (supposedly not known in 
advance). The key idea consists in splitting the original 
biphoton beam for two identical parts by a non-selective 
beam splitter (BS) and performing a series of coincidence 
photon-counting measurements in the arising channels I 
an II (see Fig. [5}i). Measurements have to be done with 
polarizers installed in the horizontal (x) and vertical (y) 
directions, as well as along the axes x' and y' turned for 
the angle 45° with respect to, correspondingly, x- and y- 
axes (Fig.[2t)). As shown below this set of measurements 
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D 

JV^ 


* 



</(V 45 ) 



«(V) 



\U- >X (H) 
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BS 



Phc. 2: (a) A scheme of and (b) two bases for coincidence 
measurements; BS denotes beam splitter, M - mirror and D - 
detectors. 



is sufficient for determining the degree of entanglement 
as well as for a full reconstruction of all parameters of 
arbitrary biphoton qutrits. 



A. Beam splitter 

In terms of biphoton wave functions, BS adds an 
additional degree of freedom to each photon - the 
propagation angles f?i and #2, which can take only two 
value, and 90°. In a nonselective BS each photon has 
equal probabilities of transmitting or being reflected. 
Thus, if the biphoton wave function in front of BS is 
after BS it takes the form 



^afterBS — * X 5(^.0 — 5g 1 90° ) {5g 2 , — ^ 2 ,90°) 



1 1 

2 I -1 



1 



(0) 



(6.1) 



As the angular variables in 'I'aftorBS are separated from 
polarization variables of \& and, besides, parts depending 
on the variables 0\ and 82 are also factorized, the angular 
factor in ^aftorBS (|6.ip does not add any additional 
entanglement to this wave function compared to "P. 
I.e., the nonselective BS itself does not increase or 
diminish the degree of entanglement of any biphoton 
states. Coincidence measurements in the state arising 
after BS are suggested here as a tool for determining 
qutrit parameters in front of BS and the degree of 
entanglement of the original biphoton state. 

As it follows from Eq. (|6.ip and from the described 
features of BS, the latter splits between the channels I 
and II only photons of a half of all incoming biphoton 
pairs, whereas another half of pairs remains unsplit, and 
these unsplit pairs are equally distributed between the 
channels I and II. Let Njj H , Nfi v , and N% v be large 
(^> 1) total amounts of biphoton pairs with coinciding 
(HH, VV) and differing (HV) polarizations of photons 
generated in a crystal per some given time. Then the total 
amount of generated pairs is N*° a t lrs = N^+Nfty+N^y 

,+„j „u„+, — „ :„ i\jtot _ 
phot 



Kftot 1 ATtot 1 ]\ftot 

and the total amount of generated photons is N* ot 



2A r '°L 5 - Among all these generated photons there are 



pairs 

ATjf = 2N% H + Ntfy horizontally and Ny 1 



N tot 



2Ny V vertically polarized photons. Relative amounts 
of horizontally and vertically polarized photons can be 
interpreted as single-particle (absolute) probabilities for 



s 



photons to have horizontal and vertical polarizations 



w 



M 
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ly phot 



1 Ajtot 

2 ly HV 



N tot 
ly phot 



N HH ' 



ATtot I ATtot 
1 ATtot 

' 2 IX HV 



(6.2) 



Though these equations are nice, unfortunately, in 
experiment amounts of generated pairs are not directly 
observable. 

After BS we get \N%£, unsplit pairs in each channel, 
and photons of the remaining ^iV^°', split pairs are 
equally distributed between the channels I and II (Fig. 
[3]). Owing to unsplit pairs results of photon counting by 



//// 



(a) 

m 



HH 



II 



II 



(c) 

w a for a photon "1" to have polarization a under the 

a' 

condition that another photon, "2", has polarization a' . 
They are defined as ratios of specific coincidence detector 
counts (16.31) to the total amounts of coincidence counts 



N, 



d(c) 



£<r,<r' N, 



d(c) 



(6.4) 



where 



Evidently, the sum of all conditional probabilities f|6.4[) 
equals unit 



= l[N%k+Ngtr+N%r] = |i\^ r . (6.5) 



(c) 

, = W H 




\-2w% ] 




(c) 
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H 
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- w v 


cr' 





1 



(6.6) 
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HV 
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Eqs. (|6.4p and (|6.5p show that, in principle, 
experimental measurement of conditional probabilities is 
straightforward because they are expressed explicitly in 



terms of the coincidence detector counts N, 



d{c) 



111 



terms of amounts of pairs in a beam expressions for 
conditional probabilities take the form 



Phc. 3: Distribution of photons in channels I and II after BS for 
(a) HH and (b) HV pairs; double lines are unsplit pairs, arcs with 
arrows indicate photons of split pairs. 

a single detector in one of the channels I or II can be 
rather confusing, because when an unsplit pair comes 
to the detector the latter can produce one click instead 
of two. In contrast, coincidence measurements by two 
detectors in channels I and II [Fig. register only 
photons of split pairs, which makes such measurements 
quite unambiguous and informative. 

Relations between amounts of single photons and the 
corresponding total amounts of pairs are different in 
the cases when single photons arise from split HH and 
VV pairs, and when they arise from split HV pairs. 
As seen from Fig. [31 the amount of H- or F-polarized 
photons arising from the split HH or VV pairs equals to 
l-^HH or h-Nyy (Fig.[3ji), whereas amounts of the same 
photons arising from the split HV pairs is equal to jNjf v 
(Fig. [3f>). From these results we easily find amounts of 

detector coincidence counts obtainable when 

a' 

a- and a' -polarizers are installed in front of detectors, 
correspondingly, in the channels I and II 



N 



d{c) 



H 



_ V ATtot 

h 2 Nhh > 



AT, 



d(c) 



_ V ATtot 

v - 2 Nyv > 



N 



d(c) 
H 



N. 



d{c) 



H 



ATtot 



(6.3) 



where r\ is the efficiency of detectors. Relative amounts of 
coincidence counts have sense of conditional probabilities 



(c) 
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N tot , x 

ATtot ' W V 
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ATtot ■ 
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2N tot 
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(6.7) 



By the definition of conditional probabilities, their sums 
can be constructed to give single-particle, unconditional 
probabilities: 



H 



w 



,(<■■) 



H 



w 



,0) 



H 



1 ATtot 

2 ly HV 



Ajtot 



ATtot 

pair 

_ 1 Ajtot 

2 Jv i?y 



ATtot 

pair 



(6.8) 



These expressions are absolutely identical to Eqs. (|6.2p 
derived simply from counting amounts of photons in 
the original biphoton beam without any beam splitter. 
Coincidence of Eqs. (|6 . 2[) and (|6.8p proves that, indeed, 
coincidence measurements after BS can be used for 
getting information about such features of the beam 
before BS as characterizing it unconditional (single- 
particle) probabilities for photons to have horizontal or 
vertical polarizations. In other words, Eqs. (|6.8|) show 
that though owing to unsplit pairs direct measurement 
of single particle probabilities after BS is problematic, 

(s) (s) 

nevertheless, the probabilities w H and w v can be found 
from results of coincidence measurements. 

Conditional probabilities of Eq. (|6.7p can be used for 
a partial reconstruction of the qutrit parameters from 
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experimental data. Indeed, directly from Eq. (|3.1j) for the 
density matrix p we find the following relations between 
the absolute values of all three coefficients Ci.2,3 in the 
wave function (|2.11[) and the conditional probabilities 



■w 



H 



\Ci\ 



w 



(c) 



|C,| 2 



w 



H 



(c) 



(6.9) 



For single-particle probabilities (|6.8[) relations with 
parameters of qutrits Cj have the form 



w 



(») 
if 



[Cil 



|C 3 | 2 



C 2 | 2 . 



(6.10) 



B. Qutrits in a turned basis 

In addition to absolute values of Cj, we have to 
find their phases, 921.2,3, and this requires additional 
measurements, e.g. with polarizers installed at 45° to 
x- and y-axes, as shown in Fig. [5f). In theory, this is 
equivalent to description of the biphoton wave function 
in a turned basis. In a general case, let x a and y a be 
axes turned for an angle a with respect to x- and y- 
axes. Let polarizations along the axes x a and y a be 
denoted as H a and V a . The corresponding one-photon 

wave functions are (q)^ and (°) " • From these one- 
photon wave functions we can construct the two-photon 
basis wave functions ^h^H", ^H a V a i an d ^vv" in the 
same way (|2.8p - (|2.10p as ^hh, ^hv, and ^vv were 
constructed from the one-photon wave functions ( J ) and 
( 1 ) . Let us present the wave function (|2.11j) in terms of 
expansion in the basis wave function of the frame (x a , y a ) 

* = Ci(a)*H«ff« +C 2 {a)^H"V +C 3 (a)^ VaVa . 

(6.11) 

One-photon wave functions in the a- and original frames 
are related with each other by the evident transformation 
formulas 



cos a 



sin a 



■ 1 \ («) 
.0) 



/0\( Q ) 

sm a ( 1 ) 
cos a (J 



.12) 



By applying these formulas to all terms and columns in 
Eq. (|2.11|) and regrouping the arising a-frame column 
products we reduce finally the biphoton wave function 
(|2.1ip to the form (|6.1ip with the following relations 
between the coefficients 61,2,3 (a) and Ci 2,3 

Ci(a) = cos 2 aC\ + \/2 cos a sin a C2 + sin 2 aC 3 , 

C2{a) = — -\/2 cos a sin a (C\ — C 3 ) + cos2aC2, 

C3 (a) = sin 2 a C\ — y/2 cos a sin a C2 + cos 2 a C3 . 

(6.13) 



Note that similar transformation formulas for the 
constants C± (|2.13j) have the form 



C+(a) = C+, 
C- (a) = cos 2a C_ + sin 2a C 2 
C 2 (a) = — sin 2a C_ + cos 2a C2 

This means, in particular, that 



(6.14) 
(6.15) 



*L Q) = cos2a*_ 



+ ! 



sin 2a4'^fv', 



= cos2avI/ H y - sin2a\l/_, 



(6.16) 



i.e. the function ^v"' has the same form in all a-frames, 

whereas the functions <I>_ and VP^y transform into each 
other with changing a. 

In a special case a = 45° Eqs. (|6. 13|) are reduced to 

Ci,3(45 ) = — ^— ± - ? =, C 2 (45 ) = — ? =-. 

(6.17) 

Similarly to (|6.9p . the qutrit parameters in the basis 
turned for 45° can be expressed in terms of the 
corresponding conditional probabilities 



(o) 



#45 



|C 1 (45°)| 2 , w% 



w 



(c) 
V 45 



H 45 



V 1: ' 



= |C 3 (45°)| 2 
^ 45 (6.18) 
i|C 2 (45°)| 2 . 

(c) 

In their turn, the conditional probabilities w L 

a <r' 45 

are related by equations identical to (|6.7p to the 
corresponding amounts of counts with polarizers installed 
in front of detectors along the directions of the x' and y' 
axes in Fig. © 



w 



(c) 



iV 



d(c) 



(6.19) 



Similarly to (|6.8p and (|6.10p we can define also single- 
particle probabilities via conditional ones in a turned 
basis 



0) 

^46 



(c) 



_ff 45 



(c) 



(a) 

V 45 



V 45 



y45 



(c) 



ff 45 



= |d(45°)| 2 , 
|C 3 (45°)| 2 . 



(6.20) 



C. Concurrence and full reconstruction of qutrit 
parameters from experimental data 

In terms of absolute values of the qutrit parameters 
61,2.3 and their phases 931.2,3, the squared concurrence C 
Eq. (|4.3p has the form 



C 2 =4|C 1 | 2 |C 3 | 2 + |C 2 | 4 
-4|C 1 ||C 3 ||C2| 2 cos(^ 1 + ^3-2^2). 



(6.21) 
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As well known and mentioned above, the common phase 
of the biphoton wave function does not affect any 
measurements. This means that multiplication of "J by an 
arbitrary phase factor e %ipo does not change any physical 
results. Let us take, for example, e lv ° = — if2, which 
makes the parameter e llpo C2 real, with corresponding 
changes of phases in two other parameters, C\ and C 3 . 
Equivalently, keeping in mind this procedure we simply 
can take C2 real and, in a general case, Ci i3 complex, 
Ci = \Ci\e 1 ^ and C 3 = |C 3 |e^ 3 . With real C 2 ((p 2 = 0), 
Eq. (|6.21|l takes the form 



|2CiC 3 - c. 



2|2 



4|C 1 | 2 |C 3 | 2 + C\ - 4|d||C 3 |C 2 cos(^ + <p a ). (6.22) 



For finding phases ipi and ip 3 we can use Eqs. (|6.1T[> . the 
second of which gives 



|C 2 (45°)p 



\G\~C 3 \ 2 



|C 1 | 2 + |C 3 | 2 



|Ci||C 3 | cos(<^i - tp 3 ). 



.23) 



One equation for two unknown phases (pi and ip 3 can be 
obtained from the difference of squared absolute values 
of the parameters Ci(45°) and C 3 (45°) in Eqs. (f6~T7f : 

|Ci(45°)| 2 - |C 3 (45°)| 2 = \/2C 2 Re(Ci + C 3 ) 
= \Z2C 2 [|Ci|cosv3i + |C 3 | costal. (6.24) 



In a general case, Eqs. (|6.23[) and (|6.24p hardly can be 
further simplified to yield a simple analytical formula 
for the squared concurrence (|6.21[) in terms of the 
experimentally measurable conditional probabilities. But 
numerical solution of Eqs. (|6.23|) and (|6.24j) with all 
parameters |C| and |C,;(45°)| known does not represent 
any problems. Thus, by getting a full set of measurements 
of coincidence counts in two bases [xy) and (x'y') (Fig. 

2b) one can find all conditional probabilities 



(ECU) and w { % r (EUD. Then, from Eqs. (JHU) and 

a' 45 

(|6.18p one finds all absolute values of parameters C 
and and Cj(45°). And, finally, with known values of |Cj| 
and I Ci (45°) I, one solves numerically Eqs. (|6.23j) and 
(|6.24|) and finds the phases <p\ and ip 3 with (f2 = 0. 
This procedure permits to reconstruct completely all 
qutrit parameters and to find its degree of entanglement 
()6.22[) from experimental coincidence measurements. 
This procedure might be referred to as an alternative 
protocol of quantum state tomography of biphoton-based 
qutrits rather than those one described in [3. Hij. 

Note, that if it's known in advance that all parameters 
d are real, the procedure of finding the degree 
of entanglement from experimental data significantly 
simplifies. In this case it's more convenient to begin 
with Eq. (|4.ip for the inverse Schmidt entanglement 
parameter, which is easily expressed in terms of single- 
particle probabilities (|6.10p . For the upper line on the 



right-hand side of Eq. (|4.1[) we get 



I Co 



C 3 I 2 + ^ 



w 



(»)2 
H 



w 



(6)2 



1 



iu H — w 



(6.25) 



On the other hand, with the help of Eqs. (|5T7)) and ([Ol , 
the lower line on the right-hand side of Eq. (|4.1[) can be 
easily expressed in terms of qutrit parameters and single- 
particle probabilities in the basis turned for 45°: 

\C* 2 C, + C 2 C* 3 \ 2 ^ CUC! + C 3 ) 2 



{[d(45°)] 2 -[C 3 (45°)] 2 }^ 



(«) 0) 



(6.26) 



Altogether this gives 



K- 1 = -{\ 



(s) 



V< 45 > 



(6.27) 



and 



C 



(6.28) 



7. QUQUARTS 
A. Definitions and wave functions 



As mentioned in the Introduction, in the generally 
accepted treatment, ququarts are considered as two-qubit 
states. In an abstract form, with unspecified pairs of 
distinguishable particles, the state vectors of ququarts 
are taken in the form [28j 

|*) 2g6 = C x |00) + C a |01) + C 3 |10) + C 4 |ll) (7.1) 

with one-qubit single-particle states |0) and |1). Written 
in the form of two-row columns, the wave function of the 
state (17.111 takes the form 



^2qb = Ci 



+c 3 



With this wave function one can easily obtain the reduced 
density matrix of Ref. [j| 

_ (\Ci\ 2 + \C 2 \ 2 c 1 c^ + c 2 ct\ 

Pr(2 qb ) - + ^ ^ + |C4|2 j , 
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Schmidt entanglement parameter 



variables 



A- 



2qb 



Tr(p 2 r2qb ) l-2|G 1 G4-C a C 3 p : 



(7.4) 



and concurrence of Ref. [2£ 



a 



2qb 



2|ClC 4 — C2C3 



(7.5) 



Though this approach and Eqs. (|7.3p ~ (|7.5p are rather 
widely used and accepted, we claim that they are 
inapplicable for biphoton ququarts formed by non- 
degenerate SPDC pairs. The first objection is the 
symmetry. The wave function (|7.2[) is asymmetric with 
respect to the permutation of particle variables 1^2, 
and this is strictly forbidden for any biphoton wave 
functions. A simple symmetrization of the wave function 
(|7.2j) is insufficient and, actually, it does not help much, 
because it reduces the wave function (|7.2[) to the qutrit's 
wave function (|2.11[) . with C2 + C3 playing the role of 
C2 in Eq. (|2.1ip and with all peculiarities of ququarts 
completely lost. 

The second objection concerns photon degrees of 
freedom and dimensionality of the Hilbert space. In 
contrast to the traditional treatment, biphoton ququarts 
formed by non-degenerate SPDC pairs of photons have a 
higher dimensionality than qutrits. Single-photon states 
from which the biphoton ququarts are constructed are 
not qubits and they form the Hilbert space of a 
dimensionality d = 4, rather than d = 2 occurring in the 
case of qutrits. Indeed, if for example, SPDC photons 
in pairs have two different frequencies w/j (high) and 
loi (low), Uh > w;, one cannot say for sure, which of 
two photons gets a higher or lower frequency. Hence, the 
photon frequency uj becomes a second photon variable, 
additional to polarization and taking two values, uj^ 
(high) and w;. In other words, each photon has now 
two degrees of freedom, polarization and frequency (3l| . 
Together, they make four states which can be occupied 
by each photon, or four combined polarization-frequency 
modes, Hh, HI, Vh, and W (instead of two polarization 
modes H and V in the case of degenerate photons). Four 
modes correspond to the dimensionality of the single- 
photon Hilbert space d = 4, in contrast to d = 2 and 
two modes in the cases of qutrits and of the traditional 
two-qubit model of ququarts. The four single-photon 
polarization-frequency state vectors are now given by 



«klO>: 4nl0>> «Ul°>. a U°) 



(7.6) 



The corresponding single-photon wave functions, 
describing states with two degrees of freedom, have two 
factors, depending on polarization (a) and frequency (u) 



8<t,H$u. 



$<T,H&u>,u>i = 



pol 



pol 



® 



pol 



pol 







(7.7) 



The superscripts po1 and u are used here for 
differentiating the polarization- and frequency- parts of 
the two-variable wave functions written in the form of 
direct products of two-row columns. 

As dimensionality of the single-photon Hilbert space is 
d = 4, these states are d=4-qudits, rather than qubits. 
Then, as well as qutrits are two-qubit states, the biphoton 
ququarts are two-qudit states. Their basis state vectors 
are given by all possible products of two different creation 
operators times the vacuum state. But, if we assume that 
two photons in each non-degenerate SPDC pair definitely 
have different frequencies, the products of two creation 
operators corresponding to coinciding frequencies have 
to be excluded to give finally only four two-photon basis 
state vectors 
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I* 
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HH I 
(qqrt), 
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O ffh°Ll >. 



a Vh a Hl 
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(qqrt)\ 
HV / 
(qqrt), 
VV I 



a Hh a Vl 
J J 



|0>, 

|0), (7.8) 



where, of course, all creation operators commute with 
each other. The basis wave functions corresponding to the 
state vectors (|7.8p are obtained with the help of general 
rules of quantum electrodynamics (|2.4p : 




(7.9) 




(7.10) 
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(7.11) 




(7.12) 



In a general form, the state vector and wave function of 
the biphoton ququart are given by superpositions of four 
basis state vectors (|7.8p and four basis wave functions 



+C 3 \^ q H rt) )+C^ q v t} ) 



(7.13) 



and 



^(qqrt) = Cl ty(W rt ) 



HH 



(7.14) 

As it should be, both all basis wave functions (|7.9[) - 
(|7.12p and the general ququart wave function (|7.14[) are 
symmetric with respect to the particle permutations 1 # 
2. Also these expressions take into account properly a 
higher dimensionality of the biphoton ququarts compared 
to qutrits and compared to the traditional simplified 
vision of biphoton ququarts based on Eqs. (j7.1[) - (|7.2D . 

Once again, dimensionality of the one-photon Hilbert 
for non-degenerate photons is d = 4. The dimensionality 
of the biphoton Hilbert space is D = d 2 = 16. 
The question is why do ququarts are characterized 
by only 4 rather than 16 basis wave functions? A 
general answer is because of some restrictions, owing 
to which 12 basis wave functions are excluded. One of 
these restrictions was mentioned above: we assume that 
frequencies of photons are always different. This excludes 
states \h, cr; h, cr') and \l, cr; I, cr'), which give rise to 6 
symmetric wave functions. Their exclusion does not have 
fundamental reasons and is related rather with the most 
often met experimental conditions. The second, deeply 
fundamental restriction is symmetry. Six remaining and 
excluded wave functions are antisymmetric, and they 
cannot be realized with photons at all. Thus, in ququarts 
all 12 excluded basis wave functions either are missing 
in all possible superpositions or can be made present 
symbolically but with obligatory zero coefficients, like 



the only antisymmetric wave function in the qutrit wave 
function (|2.16p . As well as in the case of qutrits, the 
existence of ququarts with only 4 basis wave functions 
is crucially related to the requirement of symmetry of all 
biphoton wave functions. 



B. Degree of entanglement 

As known [23| , in the case of pure two-particle states 
the dimensionality of the single-particle Hilbert space 
d determines directly the maximal achievable value 
of the Schmidt entanglement parameter, A max = d. 
Consequently, a two times higher dimensionality of 
the double-qubit Hilbert space (compared to qubits) 
doubles the maximal achievable degree of entanglement 
of ququarts if it's evaluated by the Schmidt entanglement 
parameter, -/f max = 4 (to be compared with -ftr max = 2 in 
the case of qutrits) . To find the entanglement quantifiers 
of ququarts in a general case, following to the standard 
procedure, we have to start with finding the density 
matrix p(99 rt ) , corresponding to the wave function *(«« rt ) 
(|7.14[) . Dimensionality of p^i rt > is 16 x 16 which is too 
large to be shown explicitly. But the reduced density 
matrix pf qrt ^ is much more compact, its dimensionality 
is 4 x 4 and, explicitly, it is given by 



P 



(qqrt) 



c l c 3 + c 2 C d 



/ icir + \c 2 r 

11 o |Ci| 2 + |C 3 | 

2 I cfc 3 + C%C 4 C 3 | 2 + |C 4 | 

\ cfc 2 + C% C 4 



c l c 2 + c 3 c t 


C 2 I 2 + C 4 | 2 , 



(7.15) 



The derived expression for the 4x4 reduced density 
matrix (|7.15p differs significantly from the 2x2 density 
matrix of Ref. [f| where ququarts were considered as 
two-qubit states. As explained above, the key reasons 
of this and other differences between this paper and 
Q and others are in the symmetry of wave functions 
and in a higher dimensionality of ququarts compared to 
qutrits properly taken into account in Eqs. (|7.9|) - (|7.12[) 
and (l7J4)l . 

By calculating the squared reduced density matrix 
(pr ) and its trace, we find that the Schmidt 
entanglement parameter of ququarts is simply twice 
larger than that obtained in the two-qubit model (|7.4[) 



K 



1 - 2|CiC 4 



C 2 C 3 | 2 



2K- 



2qb- 



(7.16) 



As for the concurrence, rigorously, its original Wootters' 
definition (2(| is invalid for two-qudit states with d > 2. 
But instead, one can use the so called "/-concurrence" 
[23j . by definition, determined via the Schmidt 
entanglement parameter K as C'i = a/2(1 — A' -1 ). 
Defined in this way, the /-concurrence does not provide 
any new information about entanglement of ququarts 
compared to that provided by the Schmidt entanglement 
parameter K. But the /-concurrence can be useful for 
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comparison with other entanglement quantifiers, and 
for comparison with results of the two-qubit model of 
ququarts . In a general case Eq. (|7.16| yields 



C z = ^l + ZlCiQ-CaCal 5 



(7.17) 



As it follows from Eqs. (|7.16p and (|7.17p . the maximal 
achievable values o f K and Cj for ququarts are A max = 4 
and Cj max = in agreement with the general 

expectations for qudits A max = d and Cj max = 
^2(1 - d' 1 ) [H|. The minimal values of the Schmidt 
entanglement parameter and /-concurrence are achieved 
when C1C4 = C2C3, and in this case K = K m [ D = 2 
and Cj = Cj nml = 1. This means that all biphoton 
ququarts are entangled and nonfactorable (in contrast to 
this, earlier Q , in a two-qubit model, it was assumed that 
at C1C4 = C2C3 ququarts are factorable). In particular, 
all basis states of ququarts (|7.9p - (|7.12p are minimally 
entangled. It may be interesting to note that the nature 
of entanglement of different basis states is different. 

The wave functions ^^h^ an< ^ ^vv are seen ^° ^ e 
factorized for polarization and frequency parts, and only 
their frequency parts are entangled, i.e., the states ^^h^ 
and 'fyyy^ have only a purely frequency entanglement. 

On the other hand, the wave functions $^ rt ' and 
are not factorized for polarization and frequency parts, 
and in these cases we have a nonseparable frequency- 
polarization entanglement. 

Some examples of maximally entangled ququarts are 



7f (* 



(qqrt) , j(qqrt) 

HH 1 T vv 



(qqrt) , jjqqrt) 
HV ^^VH 



Cl — ±64 — C3 — C2 — Oj , 

d = c A = 0, c 2 = ±c a = ^) , 



(qqrt) 
HH 



(qqrt) 
HV 



(qqrt) 
VH 



(qqrt) 
VV 



(7.18) 



and all functions similar to the last one but with the 
only sign "minus" located in front of other terms, , 

ftiiqrt) y(qqrt) 
* HV 01 ^VH ■ 

To see explicitly how the degree of entanglement of 
ququarts changes from K m i n = 2 to A' max = 4, let 
us consider, for example, the case C\ = cos 0, C4 = 
sin0, C3 = C2 = 0, i.e., the wave function 



4> 



sin^^. 



(7.19) 

For this state Eqs. (fTTTBj) and (|7T7|) give K(<j>) = 4/(1 + 
cos 2 20) and Cj(0) = ^l+~§ sin 2 20. These functions 
are plotted in Fig. 2] together with the subsystem entropy 
S r (</>) = — X)i=i Ajlog 2 Aj, where A 2 ; are eigenvalues of 



the reduced density matrix 
tyiiqrt) 



^4, (|7.19[) the matrix p'^^ is very simple 



|7TT5l) . In the case 



(qqrt) _ 



( \ COS 2 1 











\ COS 2 1 



^ sin 2 1 



\ sin 2 y 



(7.20) 



Evidently, its eigenvalues are Ai = A2 = | cos 2 and 



A3 = A4 = i sin 0, which yields 
S r (0) = 1 — 2 (cos 2 01og 2 I cos 0| + sin 2 01og 2 | shi( 



(7.21) 



^ ^^^v ^ 

















jt/4 



re/2 



3tc/4 



Phc. 4: The Schmidt entanglement parameter K (</)), /-concurrence 
Ci((f>), and subsystem entropy S r (4>) for the ququart characterized 
by the wave function 17. 191 . 



Eigenvalues of the reduced density matrix p^T^ are 
twice degenerate Ai = A 2 and A3 = A4. If none of 
them equals zero, the Schmidt decomposition of the 
wave function Vl/^ contains four terms of Schmidt-mode 
products. However, at = 0, when the entropy SV(0) 
is minimal, two of these eigenvalues turn zero, A3 = 
A4 = 0, and two Schmidt-mode products in the Schmidt 
decomposition disappear. It is the reason why in this 
case the degree of entanglement is minimal. However, 
the remaining two terms in the Schmidt decomposition 
indicate that in the case = the state 'J 4, is entangled 
and non-factorable. The degree of entanglement of the 
state is maximal when Ai = A 2 = A3 = A4 = 1/4, i.e., 
when all four Schmidt-mode products enter the Schmidt 
decomposition with equal weights, and this occurs at 
= 45° and = 135°. 

In a general case, ququarts are multi-parametric 
objects, and it's rather difficult to show all their features 
in a limited amount of pictures. However, it may be 
interesting to show one example more, when the ququart 
wave function has the form 



cos0^ rt) 



•sin 0^ 



(qqrt) 
HV 



(qqrt) 
VH 



V2 2 

(7.22) 

for which Ci = cos 0/V2, C 4 = sin 0/\/2, and C 2 = 
C3 = 1/2. The Schmidt entanglement parameter A(0) 
and /-concurrence C(0) of the state (|7.22[) are shown in 
Fig. [5l Asymmetry of these curves with respect to the 
substitution — > n — indicates sensitivity of the degree 
of entanglement with respect to the sign of the product 
C1C3, i.e., to phases of these real constants. 
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Phc. 5: The Schmidt entanglement parameter K ((/>), and 7- 
concurrence Cj(0) for the ququart characterized by the wave 
function ^ I j7.22^ . 



To conclude this subsection, let us describe here 
one example more of a rather peculiar polarization- 
angular ququart obtainable from qutrits with the help 
of a beam splitter and post-selection. As mentioned 
earlier, polarization- frequency and polarization-angular 
ququarts are equivalent. Nondegeneracy of photons can 
be provided by different direction of photons with equal 
frequencies in the noncollinear SPDC process. Another 
way of providing noncollinear propagation of originally 
collinear and frequency-degenerate photons is in using 
a nonselective beam splitter. As discussed above in 
Section 6A, after BS photons acquire new variables, 
9 1 and 02, which can take one of two values each, 
0° and 90°. The combined polarization-angular wave 
function after BS is given by Eq. (|6.ip . This wave 
function is not yet a complete analog of ququarts 
considered above because of unsplit pairs in Fig. [3] 
Mathematically they correspond to terms in Eq. (|6.ip 
proportional to 5g lt oSg 2 and <5ei,90° ^02,90° ■ in the case of 
frequency nondegenerate biphoton states similar terms 
would describe both photons with high or both with low 
frequencies, which we have excluded (see the discussion 
in the previous Subsection). Owing to these terms with 
coinciding angles 6\ and 02, as mentioned above, the 
angular part of the combined wave function in Eq. (|6.1[) 
does not provide any additional entanglement to the 
polarization part \& (qutrit). However, if we make a post- 
selection by killing (eliminating) in some way unsplit 
pairs in both channels I and II of Fig. [31 we get a true 
polarization-angular ququart with the wave function of 
the form 

Wirt = ^(a,a')^=[Sg u0 Sg 2 ^o+Sg ug0 oSg 2fi }, (7.23) 

where ty qtr {o,o') is an arbitrary qutrit wave function 
and, as usual, a and a' are polarization variables of 
photons. Peculiarity of the ququart, characterized by the 



wave function (|7.23|) consists in factorization of parts 
depending on polarization and angular variables. For this 
reason the total Schmidt entanglement parameter of the 
state (|7.23p appears to be factorable too: 

K qtr x K anale 



j^qqr 



(7.24) 



where K qtr and K an9le are the Schmidt entanglement 
parameters of the qutrit ^> qtr and of the angular part 
of the wave function ^ qqrt , -^[8e u o6e 2 ,m° + 8e u 9o°8e 3 ,o]- 
Actually, evidently enough, the Schmidt entanglement 
parameter of the state, characterized by the last wave 
function equals two, K anale = 2. For this reason Eq. 
([7724]) takes the simplest form: K qqrt = 2K qtr , i.e., in 
terms of the Schmidt entanglement parameter, the post- 
selection after BS doubles the degree of entanglement 
of the original qutrit. Though interesting enough, this 
last result and factorization of the Schmidt entanglement 
parameter described by Eq. (|7.24[) are specific features 
of the ququart of the form (|7.23|) . In a general case 
the degree of entanglement of ququarts, characterized 
by any entanglement quantifiers, is unseparable for 
angular/frequency and polarization parts. 



C. Measurements 

A possible way for measuring directly parameters of 
ququarts and, in particular, their degree of entanglement 
is similar to that described above for qutrits and 
illustrated by Fig. [3] the original biphoton beam has to 
be split for two channels by a nonselective beam splitter 
and a full set of coincidence measurements has to be 
done in the xy and x'y' bases of Fig. |3f). In difference 
with qutrits, in addition to polarizers one has to install 
frequency filters in front of detectors to count amounts 
of photons in given frequency-polarization Hh, HI, Vh, 
and VI. As in ququarts all basis biphoton states are of 
the type ajetj|0) with i ^ j (where i and j numerate 
above indicated modes), distribution of photons of split 
pairs in channels I and II corresponds to the type, shown 
in Fig. [3f>. For this reason, if iV*°* are total amounts 
of generated pairs with photons in modes i and j, 
the amounts of the corresponding coincidence detector 
counts are determined by equations similar to that of 
the last formula in Eq. (|6.3p : 



N 



d{c) 



N 



d(c) 



1 N tot 
4 y hJ ' 



(7.25) 



with r\ denoting the efficiency of detectors. The full set 
of data on coincidence amounts of counts can be used 
further to find conditional probabilities for a photon "1" 
to be in some mode i under the condition that the second 
photon of the same pair is in some other mode j (similarly 
to Eq. E3 ) 



,(<=) 



A, 



d(c) 



y. ■ N' 



d(c) 



(7.26) 
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(c) 



1. 



with the normalization condition Y]. ■ w, 

— 

On the other hand, the conditional probabilities are 
determined by diagonal elements of the full density 
matrix with respect to both polarization and frequency 
variables and of both photons. Directly from expressions 
(|7.9p - (|7.12p for the basis wave functions of ququarts, 
from the general expression (|7.14[) for the ququart wave 
function, and from the definition of the density matrix 
p = \J r \E r T we find the following series of relations between 
the conditional probabilities and the parameters Cj of 
ququarts 



|Ci| 



|C 2 | 2 



|C 3 | 2 
|C 4 I 2 



2»s 


HI 


= 2^] 


Hh 




VI 


= 2^1 


Hh 




Vh 


= 2u$ 


HI 


*4l 


VI 


= 2^1 


Vh 



(7.27) 



The same relations (|7.25p - (|7.27[) between amounts of 
detector counts, conditional probabilities and expansion 
coefficients of the ququart wave function can be written 
also in the x'y' basis turned for 45° with respect to 
the horizontal- vertical basis. Thus, all set of coincidence 
measurements in two bases can be used to determine all 
absolute values of ququarts parameters in both bases, 
|C| 2 and |Ci(45°)| 2 . Relations between C(45°) and Q 
are easily found in the same way as in the case of qutrits, 
and in the case of ququarts they are given by 



d(45°) = i (Cx+^ + Ca + Q), 
C 2 (45°) = | (-d + C 2 - C a + C 4 ) 
C 3 (45°) = ± (-Ci - C 2 + C 3 + C 4 ) 

C 4 (45°) = | (d - C 2 - C 3 + C 4 ) . 



(7.28) 



The squared absolute values of expressions on the left- 
arid right-hand sides of Eqs. (|7.28[) take the form of 
equations for the phases yi.2,3.4 of the constants C%. 
For shortening formulas these equations can be written 
grouped in pairs to give 



|C 1 (45°)| 2 + |C 2 (45°)| 2 -i = 
|Ci||C 3 | cos((^i - ip 3 ) + \C 2 \\C 4 \ cos(<^2 



<p 4 ), (7.29) 



|C 1 (45°)| 2 + |C 3 (45 )| 2 -i = 
IdHCal cob(v?i - <p 2 ) + |C 3 ||C 4 | cos(^ 2 - ^ 4 ), (7.30) 



|C 1 (45°)| 2 + |C 4 (45°)| 2 -i = 
+|Ci||C 4 | cosfai - ^4) + |C 2 ||C 3 | cos(^ 2 - tp 3 ). (7.31) 

Because of normalization Y^ li |C(45°)| 2 = 1, one can 
get only three independent equations of the type (|7.29[) - 
(|7.3ip from four Eqs. (|7.28[) . The fourth equation for 



finding four phases ^1,2.3,4 follows from the above 
mentioned fact that the phase of the wave wave function 
as a whole does not affect any observable 
quantities and can be arbitrarily chosen. For this reason 
one can put some additional condition restriction for 
phases, e.g., such as 



<Pl + + <P3 + <£i = 0. 



(7.32) 



Four equations (|7.29j) - (|7.32|) are sufficient for finding 
numerically all phases as soon as the squared 

absolute values |C| 2 | and |C(45°)| 2 are found from 
coincidence measurements. Thus the procedure described 
above provides a full reconstruction of the ququart 
wave function ij/(?9 rt ) (|7.14p an d> in particular, can be 
used for determining the degree of entanglement. This 
procedure serves as an alternative protocol of quantum 
state tomography of ququarts described in the work Q . 

Note that, in analogy with the case of qutrits, 
determination of entanglement quantifiers of ququarts 
from experimental data is simplified if the parameters 
Ci j2;3i 4 are known in advance to be real. In this case the 
key element of Eqs. ([7T6l) and ([7TT7P for the Schmidt 
entanglement parameter and /-concurrence, |CiC 4 — 
C2C3I 2 can be written as 

|CiC 4 — C2C3I 2 = C\C\ + C 2 C 3 — 2C1C4C2C3. (7.33) 

On the other hand, by summing the squared first and 
last lines of Eqs. (|7.28p we get 



CiC 4 + C 2 C 3 = C 2 (45°) + C 4 2 (45°) 2 



(7.34) 



Now, with a simple algebra, we find from Eqs. f|7. 33|) and 
(I7.34p the following final expression for |CiC 4 — C2C3I 2 

|dC 4 - C 2 C 3 | 2 = 2 (|d| 2 |C 4 | 2 + |C 2 | 2 |C 3 | 2 ) 
-[|C 1 (45°)| 2 + |C 4 (45°)| 2 -i] 2 , (7.35) 



where all terms on the right-hand side are expressible via 
conditional probabilities either in the original (xy) basis, 
or in the basis turned for 45° . This solves the problem of 
determining the entanglement quantifiers K or Cj . 



8. DISCUSSION 

One of the key elements of the approach of this 
paper consists in the analysis how does the obligatory 
symmetry of biphoton wave functions affect the particle 
entanglement. The conclusion is that the symmetry 
restrictions on biphoton states imposed by the Boze 
statistics are found to be important not only for correct 
calculations of the entanglement quantifiers but also for 
existence of biphoton qutrits and ququarts as classes of 
states with, correspondingly, only three and four basis 
states. Actually, for biphoton states with dimensionality 
D — 4 (qutrits) and D = 16 (ququarts) one can expect 
appearance of, correspondingly, four and sixteen basis 
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wave functions. In both cases such formations would be 
more complicated than qutrits and ququarts. To return 
to true qutrits and ququarts. we have to restrict the 
amount of basis wave functions by excluding some of 
them, and the exclusion has to result from some physics 
rather than being made artificially. The key role in such 
restriction of the amount of basis wave functions belongs 
to the symmetry requirements. In the case of qutrits the 
wave function to be excluded is the antisymmetric Bell 
state \E f_ (|2.15p . which cannot arise in any superpositions 
of purely polarization biphoton wave functions, because 
its appearance would contradict to the Boze statistics of 
photons. Alternatively, this fourth antisymmetric basis 
wave function can be included in superpositions but with 
obligatory zero coefficient, as in Eq. (|2.16p . In the case 
of ququarts the situation is even more dramatic. To get 
four basis wave functions instead of sixteen we have to 
exclude twelve wave functions. Six of them are excluded 
because we consider here only SPDC processes with non- 
coinciding frequencies of photons, and the remaining six 
basis wave functions are excluded because of symmetry 
restrictions as antisymmetric ones and, hence, forbidden 
by the Boze statistics. Thus, indeed, in the case of 
biphoton ququarts, symmetry restrictions are crucially 
important for existence of such class of states. 

The question that is often raised concerns relations 
between the fundamental entanglement arising owing 
to symmetry of wave functions and "configurational" 
entanglement arising owing to varying choice of 
coefficients C, in superpositions of basis wave functions 
(|2.1ip and (J77TJJ- As seen from our analysis, the 
entanglement arising owing to symmetry is the intrinsic 
feature of all basis wave functions of ququarts (|7.9[) - 
(|7.12p and of the qutrit's basis wave function (|2.6p . 
(|2.10p . This entanglement is fixed and cannot be changed 
by any basis transformations. On the other hand, 
entanglement arising in superpositions of basis functions 
depends on the coefficients in these superpositions. 
In this sense one can speak about configurational 
entanglement as an "addition" to entanglement arising 
owing to symmetry of wave functions. But the word 
"addition" should not be understood in this case as 
summing different types of entanglement. In a general 
case none of the entanglement quantifiers equals to a 
sum or a product of the same quantifiers, related to 
symmetry and to superposition of basis wave functions. 
I.e., the symmetry and configurational entanglement 
are unseparable from each other and together they 
determine the degree of entanglement of ququarts 
and qutrits. Related to this, it's often assumed that 
entanglement arising owing to symmetry is resourceless 
whereas the configurational entanglement is resourceful 
for applications such as protocols of quantum information 
and quantum communication like quantum teleportation 
[32l |. dense coding [33[, quantum key distribution [34| . 
etc. Moreover, sometimes this idea is extrapolated to 
the assumption that entanglement related to symmetry 
can be ignored at all or "forgotten". We think that 



such formulations are excessive. As said above, symmetry 
and configurational entanglements are unseparable from 
each other, and it's impossible to distinguish in a total 
entanglement the resourceless and resourceful parts. On 
the other hand, it's a rather interesting and important 
question what changes in applications of ququarts 
analyzed earlier in the frame of the two-qubit model when 
the latter is substituted by the above described two-qudit 
picture. We hope to return to this manifold of problems 
elsewhere. 

In connection with different opinions on entanglement 
related to symmetry of indistinguishable particles, it's 
worth mentioning the formulation by A. Peres [l8[: "We 
must now convince ourselves that this entanglement is 
not a matter of concern: No quantum prediction, referring 
to an atom located in our laboratory, is affected by the 
mere presence of similar atoms in the remote universe" . 
The second part of this statement is absolutely correct: 
far objects cannot affect any measurements performed 
solely with a laboratory object. But as for "not a matter 
of concern" , it depends on for what and under what 
conditions. For local laboratory measurements it's not 
a matter of concern, but measurements of the degree of 
entanglement assume coincidence measurements carried 
out in both laboratory and remote systems, and they 
assume a possibility of communication between systems 
for selecting only coinciding signals. Only if these 
possibilities occur, we can formulate a question about 
the degree of entanglement of the laboratory and remote 
particles. Then, if particles are indistinguishable and if 
their common state is pure, the entanglement related 
to symmetry can be a matter of concern for correct 
determining the total degree of entanglement of a 
bipartite system. To specify a little bit further these 
general speculations, let us consider the same example as 
in [l8j : two well separated identical one-electron atoms, 
"a" and "6", both in the ground states. Let the two- 
electron state be pure and let the spin state of two 
electrons be symmetric and non-entangled, Iff). Then 
the coordinate wave function of two electrons is 

i'o{xi,x 2 ) = 
[<po(x! - a)ip (x 2 -b)- ip (x 2 - a)ip (xi - b)] 

T2 ' (8,1) 

where x\ and x 2 are positions of atomic electrons, a and 
b are given positions of atomic nuclei, both x\^. and a, 
b can be vectors. We can assume that atoms are located 
far enough for the wave functions ipo (x — a) and <po {x — b) 
to not overlap, which means that J dxip* (x — a)ipo(x — 
b) = 0. Eq. (|8.1[) takes into account indistinguishability 
of electrons and the necessary Fermi-Dirac symmetry 
requirements ipo{xi,X2) = — ^0(^27 x\). Owing to this, 
we cannot say, which of electrons, "1" or "2" is located 
near the atom "a", and which near the atom "6". To 
illustrate that, indeed, local measurements with one of 
two atoms are not affected by antisymmetrization of the 
wave function, we can consider, e.g., the process of laser 
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excitation solely of an atom "a" . This situation can occur 
if the laser field is focused around the point a and is zero 
around b. The excited two-electron state with an electron 
in the atom a being at an excited level, and in the atom 
b - in the ground state has the wave function of the form: 

ip 1 (xi,x 2 ) = 
[(fi(xt - a)ip (x 2 - b)- (px(x 2 - a)Lp a {xi - b)] 

71 ' ( } 

where ip\{x — a) is the electron wave function at the 
excited level. As we don't know which electron is located 
near the atom "a" , we have to take the atom-light dipole 
interaction in the form V cosuit with V = — eeo(t)[(xi — 
a) + (x 2 — a ))]i where £o an d w are the field-strength 
amplitude and frequency of the exciting field. The rate 
of transitions found with the help of the Fermi Golden 
Rule is given by 
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(<Pi(x) I M*)) f KEi -Eq- hu), (8.3) 



where Eq and E\ are energies of the ground and excited 
atomic states. Identity of expressions in the upper and 
lower lines of Eq. (|8.3p confirms that for local excitations 
only in one of two atoms the second atom does not affect 
the rate of transitions at all and the latter can be found 
without anisymmetrization of the two-electron wave 
function, simply with (|8.ip and ipx (|8.2p substituted 

by 

Tpo = (po(x! - a)(p (x 2 - b), V>i = tpi(x! - a)tp {x 2 - b). 

(8.4) 

However, if we are interested in whether the state of 
two one-electron atoms in ground states is entangled or 
not, we have to find some of its entanglement quantifiers. 
E.g., the Schmidt entanglement parameter in the case of 
continuous variables x\ and x 2 is known [35[ to be given 
by 

A' -1 = J dxidx 2 dxidx' 2 %1>q{xx, x 2 )ipQ {x\, x 2 ) 

xi/j^xxjx'^ipofa'i,^). (8.5) 

For two remote atoms the integral on the right-hand side 
equals i, which gives K = 2 and which means that 
the state of two remote atoms characterized by a pure 
two-electron wave function is always entangled. "Always" 
means in this case as long as the state ipoixi, x 2 ) 
remains pure, i.e., as long as its purity is not hurt by 
external factors. Also this means that two atoms have 
to arise in a pure two-electron state. An example when 
it's true is the dissociation of a hydrogen molecule by 
a Fourier- limited light field for two Hydrogen atoms. 
An opposite example, when the state of electrons in 
two remote atoms is hardly pure, is when the atoms 
are obtained from different sources. In particular, the 



two-electron state hardly is pure and hardly can be 
described by a wave function when one of them is 
produced in the Earth laboratory and another somewhere 
in a remote Universe. For such atoms the quoted above 
formulation of Ref. [HI is absolutely correct: remote 
objects arising from different sources are non-entangled, 
and symmetrization or antisymmetrization of their wave 
functions has no sense because their common state 
cannot be characterized by a bipartite wave function. But 
for remote indistinguishable particles in a pure bipartite 
state their entanglement arising owing to symmetry 
can occur, can be important and measurable, and 
is unseparable from the configurational entanglement 
related to superposition of basis states. 



9. CONCLUSIONS 

To summarize, the consideration given above provided 
a systematic description of biphoton qutrits and ququarts 
and such their features as symmetry, dimensionality, 
and entanglement. As both qutrits and ququarts are 
two-bozon formations, their wave functions are obliged 
to be symmetric with respect to permutations of 
photon variables, and this condition is important for 
correct evaluation of the degree of entanglement. In our 
analysis the degree of entanglement is evaluated by such 
entanglement quantifiers as the Schmidt entanglement 
parameter, concurrence, and subsystem von Neumann 
entropy. All of them are good for pure bipartite states, 
they characterize the degree of entanglement in different 
metrics but, of course, their predictions must be and 
they are compatible with each other. The Schmidt-mode 
analysis is also applied and shown to be very fruitful 
for finding families of non-entangled, factorable qutrits. 
For qutrits we found explicitly the most general three- 
parametric families of non-entangled and maximally 
entangled states. In particular, we showed that one 
of the basis state of qutrits, in which photons have 
different polarization (horizontal and vertical), belongs 
to the family of maximally entangled qutrits. Some other 
interesting features of qutrits are analyzed, such as, 
e.g., anticorrelation of entanglement and polarization, 
and so on. In the case of biphoton ququarts, even 
more interesting, the traditional two-qubit model is 
shown to be invalid. As ququarts are produced by 
biphoton beams with nondegenerate photons (either in 
frequencies or in directions of propagation), they have 
more degrees of freedom than qutrits, e.g., polarization 
and frequency. In such case frequencies of photons are 
variables rather than simply some given numbers, and 
frequency entanglement has to be taken into account 
together with polarization entanglement. Thus, biphoton 
ququarts are shown to be two-qudit rather than two- 
qubit states, with the dimensionality of the one-photon 
Hilbert space d = 4 and dimensionality of the two- 
photon Hilbert space D — d 2 ~ 16 (in contrast with 
d = 2 and D = d 2 = 4 in the case of qutrits). This 
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new understanding of the physics of biphoton ququarts 
gave new formulas for their entanglement quantifiers 
(|7.16j) and (|7.17|) . One of qualitative consequences of 
these results is that all ququarts are entangled and 
unseparable, in contrast to earlier predictions of the two- 
qubit model. For both qutrits and ququarts, schemes of 
their complete reconstruction from experimental data are 
suggested. The schemes are based on using a nonselective 
beam splitter and carrying out full sets of coincidence 



photon-counting measurements in the usual horizontal- 
vertical basis and in the basis turned for 45°. 
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